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Abstract 

When bacteria are grown in a batch culture containing a mixture of two growth- 
limiting substrates, they exhibit a rich spectrum of substrate consumption patterns 
including diauxic growth, simultaneous consumption, and bistable growth. In pre- 
vious work, we showed that a minimal model accounting only for enzyme induction 
and dilution captures all the substrate consumption patterns (Narang, Biotech. Bio- 
eng., 59, 116, 1998; Narang, J. Theoret. Biol., accepted, 2006). In this work, we 
construct the bifurcation diagram of the minimal model, which shows the substrate 
consumption pattern at any given set of parameter values. The bifurcation diagram 
explains several general properties of mixed-substrate growth. (1) In almost all the 
cases of diauxic growth, the "preferred" substrate is the one that, by itself, supports 
a higher specific growth rate. In the literature, this property is often attributed to 
the optimality of regulatory mechanisms. Here, we show that the minimal model, 
which contains no regulation, displays the property under fairly general conditions. 
This suggests that the higher growth rate of the preferred substrate is an intrinsic 
property of the induction and dilution kinetics. (2) The model explains the pheno- 
types of various mutants containing lesions in the regions encoding for the operator, 
repressor, and peripheral enzymes. A particularly striking phenotype is the "reversal 
of the diauxie" in which the wild-type and mutant strains consume the very same 
two substrates in opposite order. This phenotype is difficult to explain in terms of 
molecular mechanisms, such as inducer exclusion or CAP activation, but it turns 
out to be a natural consequence of the model. We show furthermore that the model 
is robust. The key property of the model, namely, the competitive dynamics of the 
enzymes, is preserved even if the model is modified to account for various regulatory 
mechanisms. Finally, the model has important implications for the problem of size 
regulation in development. It suggests that protein dilution is one mechanism for 
coupling patterning and growth. 
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1 Introduction 



When microbial cells are grown in a batch culture con taining a mixt ure of 
two carbon sources, they often exhibit diauxic growth ( Monodl . 1947 ). This 



phenomenon is characterized by the appearance of two exponential growth 
phases separated by a lag phase called diauxic lag. The most well-known ex- 
ample of the diauxie is the growth of Escherichia coli on a mixture of glucose 
and lactose. Early studies by Monod showed that in this case, the two ex- 
ponenti al growth pha ses reflect the sequential consumption of glucose and 
lactose ( Monodl . Il942h . Moreover, only glucose is consumed in the first ex- 



ponential growth phase because the synthesis of the peripheral enzymes for 
lactose is somehow abolished in the presence of glucose. These enzymes in- 
clude lactose permease (which catalyzes the transport of lactose into the cell), 
/?-galactosidase (which hydrolyzes the intracellular lactose into products that 
feed into the glycolytic pathway) and lactose transacetylase (which is believed 
to metabolize toxic thiogalactosides transported by lactose permease). During 
the period of preferential growth on glucose, the peripheral enzymes for lac- 
tose are diluted to very small levels. The diauxic lag reflects the time required 
to build up these enzymes to sufficiently high levels. After the diauxic lag, 
one observes the second exponential phase corresponding to consumption of 
lactose. 

It turns out that the peripheral enzymes for lactose are synthesized only if 
lactose is present in the environment. The mechanism for the synthesis or in- 
duction of these enzymes in the prese nce of lactose and absence of glucose was 



discovered by Monod and coworkers (jJacob and Monodl . Il96ll ). It was shown 



that the genes corresponding to these enzymes are contiguous on the DNA and 
transcribed in tandem, an arrangement referred to as the lac operon (Fig. la). 
In the absence of lactose, the lac operon is not transcribed because a molecule 
called the lac repressor is bound to a specific site on the lac operon called the 
operator (Fig. lb, bottom). This prevents RNA polymerase from attaching to 
the operon and initiating transcription. In the presence of lactose, transcrip- 
tion of lac is triggered because allolactose, a product of /3-galactosidase, binds 
to the repressor, and renders it incapable of binding to the operator (Fig. lb, 
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middle). ^ 



The occurrence of the glucose-lactose diauxie suggests that transcription of lac 
is somehow repressed in the presence of glucose. Two molecular mechanisms 
have been proposed to explain this repression. 



(1) Inducer exclusion (jPostma et al.l . Il993f l: In the presence of glucose, en- 
zyme IIA^'^, a peripheral enzyme for glucose, is dephosphorylated. The 
dephosphorylated IIA^^'^ inhibits lactose uptake by binding to lactose per- 
mease. This reduces the intracellular concentration of allolactose, and 
hence, the transcription rate of the lac operon. 

Genetic evidence suggests that phosphorylated llA^^^ activates the en- 
zyme, adenylate cyclase, which catalyzes the synthesis of cyclic AMP 
(cAMP). Since the total concentration of IIA^''^ remains constant on the 
rapid time scale of its dephosphorylation, exposure of the cells to glu- 
cose causes a decrease in the level of phosphorylated IIA^'*^, and hence, 
cAMP. This reduction of the cAMP level forms the basis of yet another 
mechanism of lac r epression. 

(2) cAMP activation ( Ptashne and Gann . 20021 Chap. 1): It has been ob- 



served that RNA polymerase is not recruited to the lac operon unless 
a protein called catabolite activator protein (CAP) or cAMP receptor 
protein (CRP) is bound to a specific site on the lac operon (denoted 
"CAP site" in Fig. 1). Furthermore, CAP, by itself, has a low affinity for 
the CAP site, but when bound to cAMP, its affinity for the CAP site 
increases dramatically. The inhibition of lac transcription by glucose is 
then explained as follows. 

In the presence of lactose alone (i.e., no glucose), the cAMP level is high. 
Hence, CAP becomes cAMP-bound, attaches to the CAP site, and pro- 
motes transcription by recruiting RNA polymerase (Fig. lb, middle). 
When glucose is added to the culture, the cAMP level decreases by the 
mechanism described above. Consequently, CAP, being cAMP-free, fails 
to bind to the CAP site, and lac transcription is abolished (Fig. lb, top). 

We show below that neither one of these two mechanisms can fully explain 
the glucose-mediated repression of lac transcription. 

The following three observations contradict the cAMP activation model. 

(1) The intracellular cAMP levels during the first exponential growth phase 
(~2.5 /iM) are comparable, if not higher, than those observed during the 
second exponential growth phase (~1.25-2 /xM) (Fig. 2a). 
It follows that the repression of lac transcription in the presence of glucose 



^ A s imilar mechanism serves to induce the genes for glucose transport i Phimhridgel . 
1200.4 Fig. 4). 
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Figure 1 . Molecular mechanism f or induction and repression of the lac operon in 
E. coli ( Ptashne and Gaiml . 20021 ): (a) Structure of the lac operon. The lacZ, lacY, 



and lacA genes code for /9-galactosidase, lactose permease, and lactose transacety- 
lase, respectively. The operator, promoter, and CAP site denote the DNA sequences 
which bind the repressor, RNA polymerase, and CAP-cAMP, respectively, (b) The 
states of the lac operon in the presence of glucose or/and lactose. The repressor and 
CAP-cAMP complex are denoted rep and CAP, respectively. 
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Figure 2. Repression of lac transcription in the presence of glucose is not due to 



reduced cAMP levels I Inada et al 



Kimata et al 



1997n . OD denotes optical 
density, and Miller units are a measure of /3-galactosidase activity, (a) Growth of the 
wild-type strain, E. coli W3110, on glucose + lactose. The intracellular cAMP levels 
are comparable during the two exponential growth phases, (b) Growth of E. coli 
W3110 on glucose + lactose in the presence of 5 mM cAMP. Despite the high cAMP 
concentration, /3-galactosidase synthesis is repressed during the first exponential 
growth phase, (c, d) Growth of E. coli Acya crp* and PR166 on glucose + lactose. 
The lac transcription rate in these strains is independent of the cAMP level, but 
/3-galactosidase synthesis is repressed during the first exponential growth phase. 

is not due to lower cAMP levels. ^ 
(2) When the culture is exposed to large concentrations (5 mM) of exoge- 
nous cAMP, the diauxic lag vanishes, but the lac operon still fails to be 
transcribed during the first exponential grovv^th phase (Fig. 2b). 



Excess cAMP fails to relieve the repression of transcript ion during growth o f E. coli 



on other pairs of substrates, such as glucose + melibiose ()Oka,da et al 
and glucose + galactose (see Fig. 9a of this work). 



mil. Fig. 4) 
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The disappearance of the diauxic lag implies that an elevated level of 
intracellular cAMP does stimulate the lac transcription rate. However, 
it fails to relieve the repression of lac transcription in the presence of 
glucose. 

(3) Diauxic growth persists in cells which transcribe the lac operon at a rate 
that is independent of cAMP levels. This has been demonstrated with 
two types of cells (Fig. 2c, d). In E. coli Acya crp* mutants, crp, the gene 
coding for CAP, is mutated such that CAP binds to the CAP site even 
in the absence of cAMP. In E. coli PR166, the lac promoter is mutated 
such that RNA polymerase binds to the promoter even if there is no CAP- 
cAMP at the CAP site. In both cases, transcription of lac is independent 
of cAMP levels. Yet, /3-galactosidase synthesis is still repressed during 
the first exponential growth phase. 



These results show that higher cAMP levels do stimulate the lac transcrip- 
tion rate. Indeed, the 5-fold increase in cAMP levels at the end of the first 
exponential growth phase in Fig. 2a is characteristic of cells expo sed to low 



concentrations (0.3 mM) of glucose (jNotlev-McRobb et al.L I1997I ). and it is 
likely that this serves to reduce the length of the diauxic lag. However, lac 
transcription is repressed in the presence of glucose even if the ability of cAMP 
to influence lac transcription is abolished. 



The persistence of the glucose-lactose diauxie in cAMP-independent cells has 
led to the hvpothe sis that inducer exclusion alone is respo nsible for inhibiting 



lac transcription (|Ina.da et al. I. Il 994 iKimata et all . 119971 ) . However, inducer 



exclusion exerts a relatively mild effect on lactose uptake. In E. coli ML30, 
the activity of lactose perm ease is inhibited no more th an ~40% at saturating 
concentrations of glucose ()Cohn and Horibata . h 959l Table 2). This partial 
inhibition by inducer exclusion cannot explain the almost complete inhibition 
of lac transcription. 



Thus, despite several decades of research, no molecular mechanism has been 
found to fully explain the glucose-lactose diauxie in E. coli. 



In the meantime, microbial physiologists have accumulated a vast body of work 
showing that diauxic growth is ubiquitous. It has been observed in diverse mi- 
crobial species on many pairs of substitutable substrates (i.e., substrates that 
satisfy the same nutrient requireme r its) including pairs of carbon sou rces (|Egli, 
199.4lHyder and Diikhu izenl. Il982 ; jKovijirova-Kovar and Eglil . Il998|l. nitrogen 

hosphorus sources (jPaughton et al 



sourc es (INeidhardt and Mag a saoik 
Il979h . and electron acceptors 



Pt all 



These studies show that there 
is no correlation between the chemical identity of a compound and its ability 
to act as the preferred substrate. For instance, during growth on a mixture 
of glucose and an organic acid, enteric bacteria, such as E. coli, prefer glu- 
cose, whereas soil bacteria, such as Pseudomonas and Arthrobacter, prefer the 
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(a) 



(b) 



Figure 3. Simultaneous consumption of substrates in batch cultures: (a) Simulta- 
neous consumption of fumarate (FUM) and pyruvate (PYR) during batch growth 
of E. coli K12 (c denotes the cell density in gms dry weight per liter). The sin- 
gle-substrate maximum specific growth rates on fumarate and pyruvate are 0.41 h~^ 
and 0.28 h~^ respectivelv. This gr owth pattern is observed with several pairs of 
organic acids i Narang et al. . 1997bh . (b) Simultaneous co nsumption of gluco se and 
mannitol (MTL) during batch growth of E. coli strain 158 (|Lengeler and Li 
There is significant uptake of mannitol during the first 4 hours even though the cells 
are precultured on glucose. 

organic acid ( Harder and Dijkhuizen . 19761 19821 ). However, there is a correla- 
tion between the maximum specific growth rate on a compound and its ability 
to act as a preferred substrate. 



In most cases, although not invariably, the presence of a substrate per- 
mitting a higher growth ra te prevents the utilization of a second, 'poorer', 
substrate in batch culture (jHarder and Dijkhuizen . ll982L p. 461). 



This remarkable correlation, which is reminiscent of anthropomorphic choice, 
is often rationalized by appealing to teleological (design-oriented) arguments. 
Harder & Dijkhuizen assert, for instance, that consumption of lactose is abol- 
ished in the presence of glucose because this prevents "unnecessary synthesis 
of catabolic enzymes in cells that alr eady have available a carbon and en- 
ergy source that allows fast growth" (iHarder and Dijkhuizen . ll982L p. 463). 
However, there is no mechanistic explanation for this correlation. 



Although the diauxie dominates the literature on mixed-substrate growth, 
there is ample evidence of nondiauxic growth. In E. coli K12, seve ral pairs 



of organic acids are consumed simultaneously ( Narang et al. . 1997bh . 



one ex- 



ample of which is shown in Fig. 3a. The maximum specific growth rates on 
these organic acids are in the range 0.28-0.44 h~^, which are low compared to 
the largest maximum specific growth rate sustained in a minimal (synthetic) 
medium (0.74 h~^ on glucose). Similar behavior has been observed in other 
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species, leading Egli to conclude that 



Especially combinations of substrates that support m edium or low maxi- 
mum specific growth rates are utilized simultaneously (Edi . E995i . p. 325). 

However, a closer look at data suggests that low or medium growth rates are 
not necessary for simultaneous consumption. This is evident from Monod's 
early studies with the so-calle d "A-su g ars," namelv i glucose, fructose, man- 
nitoL mannose, and sucrose (iMonodl . I1942L Il947f l.^ He found that in E. 
coli and B. subtilis, these sugars supported comparable maximum specific 
growth rates, but there was no diauxic lag during growth on a mixture of 
glucose and any one of the other A-sugars. Subsequent studies have confirmed 
that in some of these cases, both the substrates are consumed simultane- 
ously (Fig. 3b). Now, in all the cases of simultaneous consumption described 
above, the single-substrate growth rates were comparable. Thus, it is conceiv- 
able that simultaneous consumption occurs whenever the ratio of the single- 
substrate growth rates is close to 1. It turns out that this condition may 
be necessary, but it is certainly not sufficient. Although the growth rates of 
Propionibacterium shermanii on glucose and lactate are id entical fO-lfl ^^. ^ 
0.142 h~\ respectively), lactate is consumed preferentially ()Lee et al.L Il974r i. 
Similarly, the growth rates of E. coli ML308 on glucose and fructose are com- 
parabl e (0.91 and 0.73 h~^, re spectivelv). but glucose is consumed preferen- 
tially ( Clark and Holmj . 197fil ). ^ Thus, current evidence suggests that the ex- 
istence of comparable single-substrate growth rates is, perhaps, necessary, but 
not sufficient, for simultaneous consumption. It seems desirable to understand 
the mechanistic basis of this observation. 



In addition to simultaneous substrate utilization, there is some evidence that 
the substrate utilization pattern can depend on the history of the precul- 
ture. Hamilton & Dawes were among the first to observe such behavior dur- 
ing t he growth of Pseudomonas aeruQinosa on a mixture of citrate and glu- 
cose ( Hamilton and Pawed . ll959L Il96fll . 1961 ). Cells precultured on citrate 
showed diauxic growth with citrate as the preferred substrate, whereas cells 
precultured on glucose consumed both citrate and glucose. We observed a simi- 
lar substr ate consumption patte rn during growth of E. coli K12 on glucose and 
pyruvate ( Narang et al. . 1997bl ). An entirely different preculture-dependent 
pattern was obtained during the growth of a pseudomonad on glucose and 
phenol (iPanikov . I1 99,4 Chap. 3 p. 181). When the cells were precultured on 



^ It was found later that all the A-sugars are transported by the phosphotransferase 
system (PTS) dRosemaTi and MeaHoJ Fflflfl). 



^ The absence of the diauxic lag, observed in Monod's earlier studies with glucose 
fruct ose mixtures, is due to r apid de novo synthesis of the PTS enzymes for fruc 



tose (|Clark and Holmsl . Il976l . Figs. 4-5). Thus, preferential consumption without 
a lag does not imply the existence of new molecular mechanisms — it can be a 
consequence of rapid induction kinetics. 
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glucose, there was preferential consumption of glucose. Immediately after the 
exhaustion of phenol, when the cells were fully adapted to phenol, the medium 
was supplemented with additional glucose and phenol. Once again, there was 
diauxic growth, but phenol, rather than glucose, was the preferred substrate. 
In earlier work, we have argued that preculture-dependent growth patterns 
may be quite common — the lack of such data reflec ts the fact that the ef fect 
of preculturing was not investigated in most studies (jNarang et al. i ll997bl l. In 
order to facilitate their identification, it seems appropriate to determine the 
feasible preculture-dependent growth patterns. 



The goal of this work is to seek mechanistic answers for the following questions 

(1) In diauxic growth, why is the maximum specific growth rate on the pre- 
ferred substrate higher than that on the less preferred substrate? 

(2) Under what conditions are the substrates consumed simultaneously? 

(3) What types of preculture-dependent growth patterns are feasible? 

There are numerous mechanistic models of mixed-substrate growth. Many of 



lactose diauxie in E. coli JKremliuff et al.. 2001: 


ISantillan and Mackev. 2004: 


van Dedem and Moo-Yound. Il97rl Wong et al.l. 


1997 


). These models cannot 



address the above questions, which are concerned with the general proper- 
ties of mixed-substrate growth. Thus, one led to consider the more general 
models accounting for onlv those processes tha t are common t o most svstems 
of mixed-substrate growth JBrandt et al. . 2004 : Narang et al. . 1997a : Naxang . 



1 998al : iThattai and ShraimanL I2 



. Recently, we have shown that these gen- 
eral models are s imilar inasmuc h as the enzymes follow competitive dynamics 
in all the cases ( Narang . 2006[ ). However, the model in Brandt et al cannot 



capture nondiauxic growth, and the model in Thattai & Shraiman treats the 
specific growth rate as a fixed (constant) parameter, an assumption that is 
not appropriate for describing the growth of batch cultures. 



In this w ork, we address t he questions posed above by appealing to the minimal 
model in lNarane) (jl998al ). This model accounts for only enzyme induction and 
dilution, the two processes that occur in almost all systems of mixed-substrate 
growth. Yet, it captures all the batch growth patterns described above, and 
its extension t o continuous cu ltures predicts all the steady states observed 
in chemostats ( Narang . 1998bl ). Here, we show that the minimal model also 
provides mechanistic explanations for the foregoing questions. Specifically, we 
find that 



(1) If the induction kinetics are hyperbolic, the maximum specific growth 
rate on the preferred substrate is always higher than that on the less 
preferred substrate. The manifestation of this correlation in a minimal 
model containing no regulation suggests that its existence is not due 
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CO 




CO 



Figure 4. Kinetic scheme of the minimal model ( Narangl . 1998ah . 



to goal-oriented regulatory mechanisms, an ass umption that lies at the 



heart of models based on optimality principles (iMahadevan et all [2 



erty resulting from the kinetics of enzyme induction and dilution. We also 
find that the correlation can be violated when the induction kinetics are 
sigmoidal, and that the dynamics of these offending cases are consistent 
with the data in the literature. 
(2) The existence of comparable single-substrate growth rates is not suffi- 
cient for simultaneous consumption. This agrees with the data described 
above. However, we find that this condition is not necessary either. This 
is because the occurrence of simultaneous consumption depends not only 
on the relative growth rates, but also on the saturation constants for 
induction. If these saturation constants are small, there is simultaneous 
consumption, regardless of the relative growth rates. 

We show, furthermore, that the classification of the substrate consumption 
patterns predicted by the model explains the phenotypes of several mutants. 
The most striking phenotype is the reversal of the diauxie, wherein both the 
wild-type and the mutant strains display diauxie growth, but consume the 
substrates in opposite order. This phenotype cannot be explained in terms of 
the standard molecular mechanisms, but turns out to be a natural consequence 
of the minimal model. 



2 The model 



Fig. 4 shows the kinetic scheme of the minimal model. Here, Si denotes the 
i^^ exogenous substrate, Ei denotes the transport enzyme for Si, Xi denotes 
internalized Si, and C~ denotes all intracellular components except Ei and 
Xi (thus, it includes precursors, free amino acids, and macromolecules). 
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In this work, attention will be confined to growth in batch cultures. We assume 
that 



(1) The concentrations of the intracellular components, denoted Cj, Xj, and 
c~, are based on the dry weight of the cells (g per g dry weight of cells, 
i.e., g gdw~^). The concentrations of the exogenous substrate and cells, 
denoted Si and c, are based on the volume of the reactor (g/L and gdw/L, 
respectively). The rates of all the processes are based on the dry weight of 
the cells (g gdw~^ We shall use the term specific rate to emphasize 
this point. 

The choice of these units implies that if the concentration of any intra- 
cellular component, Z, is 2; g gdw~^, then the evolution of z in batch 
cultures is given by 

dz ^ _ fl dc\ 
dt ^ \c dt J 



where r+ and r~ denote the specific rates of synthesis and degradation 
of Z in g gdw~^ h~^. 

(2) The transport and peripheral catabolism of Si is catalyzed by the "lumped" 
system of peripheral enzymes, Ei. The specific uptake rate of Si, denoted 
rs,i, follows the modified Michaelis-Menten kinetics, r^^j = Vs,ieiSi/{Ks,i + 

Si). 

(3) Part of the internalized substrate, denoted Xj, is converted to C~ . The 
remainder is oxidized to CO2 in order to generate energy. 

(a) The conversion of Xj to C~ and CO2 follows first-order kinetics, i.e., 

(b) The fraction of Xj converted to C~ is a constant (parameter), denoted 
Yi. Thus, the specific rate of synthesis of C~ from Xj is Fjr^^^j.^ 

(4) The internalized substrate also induces the synthesis of Ei. 

(a) The specific synthesis rate of Ei follows Hill kinetics, i.e., Ve^i = 
Ve,ix"'/(-^e^i + where = 1 or 2. Kinetic analysis of the data 
shows th at enzvme induction ca n be hyperbolic (rij = 1) or sigmoidal 
{ui = 2) (IVagil and Yagillll 97ll l. 

By appealing to a molecular model of induction, we can express rii, 
Ve,i, and Ke^i in terms of the parameters associated with repressor- 
operator and repressor-inducer binding. It is shown in Appendix A 
that the Yagil & Yagil model of induction implies that is the 
number of inducer molecules that bind to 1 repressor molecule. Fur- 



The so-called conservative substrates, such as nitrogen and phosphorus sources, 
are completely assimilated (as opposed to carbon sources, which are partially oxi- 
dized to generate energy). During growth on mixtures of such substrates, 1^ = 1 for 
both the substrates. 
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thermore, if the enzyme is inducible, 

Kxi 

Ve,i = l^e,iOt,i, = -J^rt,i, (1) 

where Ve,i is the enzyme synthesis rate per unit mass of operator; 
Ot,i,rt^i are the total concentrations of the operator and repressor (g 
gdw^^), respectively; and K^^i, Ko,i are the dissociation constants for 
repressor-inducer and repressor-operator binding, respectively. 

(b) The synthesis of the enzymes occurs at the expense of the biosyn- 
thetic constituents, C~ . 

(c) Enzyme degradation is negligibly small. 

Given these assumptions, the mass balances yield the equations 

H " ''''' = '''''' K~+7i' ^ ^ 

dxi 1 1 dc\ _ 

dci 1 1 dc\ _ x"-* 

Itt " ""^'^ ~\cJt) ''''' = ^'XH^^"^' 

^ = (r.r., + y.r.,)-K. + re,)-(igc-. (5) 

It is shown in Appendix B that since Xi + ^2 + ei + 62 + c~ = 1, Eqs. (3)-(5) 
implicitly define the specific growth rate, denoted r^, and the evolution of the 
cell density via the relations 

dc ^ 2 

— = TgC, Tg = Y,rs,i - 1^(1 - Yi}r^^i. (6) 

i=i 1=1 

Furthermore, since Xi is small, it attains quasisteady state on a time scale of 
seconds, thus resulting in the simplified equations 

ds- 

Xi ~ , {Q} 

dt - "^'^ "^'^'^ "^^'^ ~ "^^'k:^ + [eisJ{K,, + 5,)]"« ' - ''^■V,,' 

dc 

— = TgC, Tg ^ Yirs,i + Y2rs,2, (10) 
c~ = 1 - xi - X2 - ei - 62, (11) 

where (8) is obtained from the quasisteady state relation, i.e., ~ rs,i — rx,i- 
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We are particularly interested in the dynamics of the peripheral enzymes dur- 
ing the first exponential growth phase, since it is these finite-time dynamics 
that determine the substrate utilization pattern. If the peripheral enzymes for 
one of the substrates vanish during this period, there is diauxic growth; if the 
peripheral enzymes for both substrates persist, there is simultaneous substrate 
utilization. 

It turns out that the motion of the enzymes during the first exponential growth 
phase is governed by only two equations. To see this, observe that during the 
first exponential growth phase, both substrates are in excess, i.e., Sj ^ Kg i. 
Hence, even though the exogenous substrate concentrations are decreasing, the 
transport enzymes see a quasiconstant environment (si/{Ks^i + Si) and 
approach the quasisteady state levels corresponding to exponential (balanced) 
growth. This motion is approximated by the equations 



dt "'^KZ + el 

n 

2 62 

-"2 

e,2 "T ^2 



^ = V,,2 ^nfl^n, - iYlVs,ie^ + Y2Vs,2e2) 62, (13) 



obtained from (9) by replacing Si/{Ks^i + Sj) with 1. We shall refer to these 
as the reduced equations. It should be emphasized that the steady states of 
the reduced equations ar e quasisteady states of the full system of equations 



(see iNarang et al.L Il997al for a rigorous derivation of the reduced equations) . 



The reduced equations are formally similar to the equations of the standard 
Lotka-Volterra model for two competing species, namely, 

^ = riiVi(l - aniVi - aiaiVa), (14) 

^ = r2A^2(l - a2iA^i - a22A^2), (15) 

where Ni is the population density of the i^^ species, is the (unrestricted) 
specific growth rate of the i*^ species in the absence of any competition, and 
aii,ai2 are parameters that quantify the reduction of the unrestricted spe- 



cific growth rate due to intra- and inter-specific competition ()Murravl . uS 
Thus, enzyme induction is the correlate of unrestricted growth, and the two 
dilution terms are the correlates of intra- and inter-specific competition. In 
what follows, we shall constantly appeal to this dynamical analogy. 

The dynamics of the standard Lotka-Volterra model are well understood. In- 
deed, the bifurcation diagram of the model is completely determined by the 
two dimensionless parameters, 621 = «2i/oii and 612 = 012/022 (Fig. 5). These 
parameters characterize the extent to which each species inhibits the other 
species relative to the extent to which it inhibits itself. Both species coex- 
ist precisely when they inhibit themselves more than they inhibit the other 
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Figure 5. Classification of tlie global dynamics for the standard Lotka-Vblterra 
model. The full and open circles show stable and unstable steady states, respec- 
tively. 

species, i.e., 621,^12 < 1- Under all other conditions, coexistence is impossi- 
ble. If the interaction between the species is asymmetric (621 < l,&i2 > 1 
or 621 > 1,^12 < 1), one of them is rendered extinct (species 1 and 2, re- 
spectively). If both species inhibit the other species more than they inhibit 
themselves, i.e. &2i,&i2 > 1, the outcome depends on the initial population 
densities. 

Given the formal similarity of the reduced equations to the Lotka-Volterra 
model, we expect them to display "extinction" and "coexistence" dynamics. 
Importantly, these dynamics have simple biological interpretations. Extinction 
of one of the enzymes corresponds to diauxic growth, and coexistence of both 
enzymes corresponds to simultaneous consumption. It is therefore clear that 
the bifurcation diagram for the reduced equations is a useful analytical tool. 
It furnishes a classification of the substrate consumption patterns, which can 
then be used to systematically address the questions posed in the Introduction. 
Our first goal is to construct this bifurcation diagram. 

To minimize the number of parameters in the bifurcation diagram, we rescale 
the reduced equations by defining the dimensionless variables 

Bi I 

Q = , ^=, T = t\ Ve lYiVs 1 . 

The choice of the reference variables in this scaling is suggested by the following 
fact: \JVe^i/ (V'jKj^j) and ^K^^l^V^ are upper bounds for the enzyme level and 
maximum specific growth rate attained during single-substrate exponential 
growth on saturating concentrations oi Si. Indeed, under these conditions, the 
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mass balance for Ej becomes 



= V YV 



Hence, Cj < and the maximum specific growth rate on Si, de- 

noted r™f'', satisfies the relation 



The above scaling yields the dimensionless reduced equations 

|^ = ^^T^-(ei + a6,)e„ (16) 



de2 €2 



"2 



with dimensionless parameters 



Q «, , no - (ei + ae2)e2, (17) 



'^j — — / — J^e,ii^x,i\l 77—77— = — \i 77—77—, (^ioj 



a= \ 19 



These dimensionless parameters have simple biological interpretations. We 
can view Kj as a dimensionless saturation constant for induction, and a as a 
measure of the maximum specific growth rate on S2 relative to that on 5*1. 



3 Results 



We wish to construct the bifurcation diagram for Eqs. (16)-(17). Si nce limit cv- 

cles a re impossible in Lotka-Volterra models for competing species ( Hirsch and Smal5 



19741 ) ■ it suffices to determine the steady states and their stability. 



Eqs. (16)-(17) admit at most four types of steady states: The trivial steady 
state (ei = £2 = 0), the semitrivial steady states (ei > 0, €2 = and ei = 
0,62 > 0), and the nontrivial steady state, ei,e2 > 0. We denote these steady 
states by £'00, -E107 -^01 1 and En, respectively. 

We shall consider two cases: rii = n2 = 1 and rii = 2, n2 = 1. The second case 
will serve to show the qualitative changes engendered by sigmoidal induction 
kinetics. 
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(a) 



(b) 



Figure 6. Bifurcation diagrams for the case ni = 77-2 = 1 at (a) fixed K2 > and 
(b) fixed Ki > 0. In the phase portraits, the nullchnes for ei and £2 are represented by 
full and dashed lines, respectively; stable and unstable steady states are represented 
by full and open circles, respectively. The graphs of a*, a*, and Og are represented 
by blue, red, and dashed brown curves. 



3. 1 Case 1 (rii = n2 = I) 



In this case, the scaled equations are 



dei 

H 

de2 
'dt 



Ki + ei 

£2 



a- 



K2 + £2 



- (ei + ae2)ei, 
(ei + ae2)e2. 



The bifurcation diagrams for these equations are shown in Fig. 6. They were 
inferred from the following facts derived in Appendix C. 

(1) The trivial steady, Eqq, always exists (for all a, ki,K2 > 0), but it is 
always unstable (as a node). 

(2) The semitrivial steady state, Eiq, always exists. It is (uniquely) given by 



ei = 7. , £2 = 0, 



and is stable (as a node) precisely if eil^,^^ exceeds a/K2, the ei-intercept 
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of the nontrivial nullcline for €2- ^ That is 



ei\Eio> — ^ Oi < a^{Ki,K2) = ^ O ~- 

(3) The steady state, -Eqi, always exists. It is given by 



-K2 + J/ti + 4 

ei = 0, 62 = ^ , 

and it is stable (as a node) precisely if exceeds 1/(q;ki), the €2- 

intercept of the nontrivial nullcline for ei, i.e., 

1 2 

^2^0 > ^ a > a*(Ki,K2) = 7 X ■ (21) 

(4) The surface of a*(Ki, ^2) lies below the surface of q*{k\, K2), i.e., 

a^{Ki,K2) < a*{Ki,K2) (22) 

for all Ki, K2 > 0. The notation was chosen to reflect this fact: The func- 
tions, a^{K,i,K,2) and q;*(ki,K2), represent the lower and upper surfaces 
of the bifurcation diagram. 

(5) The steady state, i^n, exists if and only if both Eiq and Eqi are unstable, 
i.e., 

a*(/ti, K2) < a < a*{Ki, K2). (23) 
It is unique and stable whenever it exists. 

The bifurcation diagrams imply the following classification of the substrate 
utilization patterns. 

(1) If a < a;*(Ki,K2), only Eiq is stable, which corresponds to preferential 
consumption of Si. 

(2) If a*(Ki, K2) < < a*(Ki, K2), only En is stable, and there is simultane- 
ous consumption of Si and S'2. 

(3) If a > «*(/€!, K2), only £^01 is stable, which corresponds to preferential 
consumption of S'2. 

Thus, the surfaces of a*(Ki, ^2) and a*(Ki, K2) delineate the boundaries of the 
substrate consumption patterns. ^ 



^ The nuUdines for refer to the locus of points on the eie2-plane at which dei/dt = 
0. In the case of the reduced equations, the nullclines for ej consist of two curves. One 
of these curves is the trivial nullcline, ej = 0; the other curve is called the nontrivial 
nullcline. 

® An analogo us classification i s also obtained when the model is extended to contin- 
uous cultures ( Naxa.ng . 19!j8hl . Fig. 10). However, the control parameters consist of 
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(a) (b) 

Figure 7. The bifurcation diagram obtained when a is increased at ki = K2 = 1. The 
red, blue, and black curves represent the locus of i?io, -Eoi, and E'n, respectively. 
The curves are full (resp., dashed) if the steady state is stable (resp., unstable). The 
bifurcation points at a = a*(l,l) = 0.62 and a = a*(l,l) = 1.62 are represented 
by full circles. 



If the point, (Ki,ft2,a), crosses either one of these boundaries, there is an 
abrupt transition in the substrate consumption pattern due to transcritical 
bifurcations. This becomes evident if a is increased at any fixed ki, K2 > 
(Fig. 7). At a = a*(/€i,K2), the substrate consumption pattern switches from 
preferential consumption of Si to simultaneous consumption of Si and 5*2 
through a transcritical bifurcation in v^hich Eiq (red curve) yields its stability 
to Ell (black curve). As a is increased further, there is another transition 
at a = a*{Ki,K2) wherein simultaneous consumption switches to preferen- 
tial consumption of S2 via a transcritical bifurcation involving the transfer of 
stability from En (black curve) to Eqi (blue curve). 

We gain intuitive insight into the bifurcation diagram by considering two lim- 
iting cases. Fig. 6 shows that if ki or K2 are large, the curves for a* and a* 
converge, and simultaneous consumption is virtually impossible. In contrast, 
if both Ki and K2 are small, there is simultaneous consumption for almost 
all a. To understand these limiting cases, observe that when ki, K2 are large, 
Eqs. (16)-(17) are approximated by the equations 

— ^ — ei (1 - KiCi - aKie2) ei, 
at Ki 

de2 « / K2 \ 

— ^ —£2 1 ei - K2e2 €2, 

at K,2 \ a J 



the dilution rate and feed concentrations (rather than the physiological parameters, 
a, Ki, and K2). 
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which are formally identical to the standard Lotka-Volterra model with an = 
^1, = aKi, = K,2/ct, and 022 = k,2- However, there is an important 
difference. The parameters, 621 = 021/aii, 612 = 012/022, are not independent 
since 621 = H2/{o^i^i) = 1/&12- But if 612 and 621 are restricted to the curve 
&21&12 = 1, Fig. 5 implies that coexistence (i.e., simultaneous consumption) is 
impossible: Ei becomes extinct if 621 < 1, and E2 becomes extinct if 621 > 1- 
On the other hand, if ki, K2 are small, the enzyme synthesis rate is essentially 
constant (quasi-constitutive). The enzymes therefore resist extinction, and 
coexist for almost all a. 



3.1.1 Dependence of substrate consumption pattern on genotype 

In the experimental literature, the influence of the physiological parameters is 
often studied by altering the genetic make-up (genotype) of the cells, and ob- 
serving the resultant change in the substrate consumption pattern (phenotype) 
of the cells. We show below that the bifurcation diagrams are consistent with 
the phenotypic changes observed in response to various genotypic alterations. 

Before doing so, however, it is useful to note that in all the experiments de- 
scribed below, the phenotype of the wild-type strain is preferential consump- 
tion of a substrate (glucose, in most cases). Since Eqs. (16) and (17) are 
formally the same, there is no loss of generality in assuming that the preferred 
substrate is 5*2, and the parameters, /€i,/€2,a, for the wild-type strain lie in 
the region, a > a* (above the red curve in Fig. 6). 

We begin by considering the cases in which the genetic perturbation trans- 
forms the substrate consumption pattern from preferential to simultaneous 
consumption. 

In wild-type E. coli, transcription of lac is abolished in the presence of glucose. 
However, mutants with lesions in the lac operator svnth esize /3-galactosidase 
even in the presence of glucose ( Jacob and Monodl . 1961 ). Thus, the mutation 



transforms the substrate consumption pattern from preferential consumption 
of glucose to simultaneous consumption of glucose and lactose. The very same 
phenotypic change is also observed in mutants with a defective lad, the gene 
encoding the lac repressor ( .Tacob and Monodl . 1961 ). To explain these pheno- 



typic changes in terms of the model, observe that mutations in the lac operator 
or lad impair the lac repressor-operator binding, i.e., they increase the dis- 
sociation constant, Ko^i. It follows from Eqs. (18)-(19) that ni decreases at 
flxed K2 and a. Inspection of Fig. 6a shows that such a change can shift the 
substrate consumption pattern from preferential consumption of 5*2 to simul- 
taneous consumption. 

If lacY, the gene encoding lactose permease, is overexpressed in E. coli PR166, 
synthesis of /3-galactosidase persists in the presence of glucose (Fig. 8a). Now, 
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Figure 8. In E. coli PR166, synthesis of /3-galactosidase persists in the presence 
of gluco se if (a) lacY is pverexpressed or (b) err, the gene for enzyme II A^'*^, is 
deleted (jKimata et al. . 1997i ). The mutant cells in (b) grow on glucose despite the 



absence of IIA^'*^ because of slow uptake of glucose by the PTS enzymes for mannose. 

in the model, overexpression of lacY corresponds to higher Ki. It follows 
from Eqs. (18)-(19) that Ki,a decrease at fixed k,2, and Fig. 6a implies that 
the observed phenotype is indeed feasible. 

In E. coli PR166, /3-galactosidase is synthesized despite the presence of glu- 
cose if err, the gene for enzyme II A^'*^, is deleted (Fig. 8a). Similarly, in the 
wild-type strain, E. coli K12 W3110, glucose is consumed before galactose. 
However, mutants with lesions in a gene encoding a transport enzyme for glu- 



cose consume the two substrates simultaneously (jKamogawa and Kurahashi 



19671 ). In these cases, the effect of the mutation is to decrease V"s,2, so that K2 
and a decrease at fixed ki. It follows from Fig. 6b that such a change could 
lead to simultaneous consumption of the substrates. 

Now, all the mutant phenotypes discussed above can be explained just as 
well by alternative hypotheses appealing only to the molecular mechanisms. 
Indeed, the first case is obviously due to impaired repressor-operator binding, 
and one can argue that the remaining two cases are due to diminished inducer 
exclusion. However, the next two examples, which involve the reversal of the 
diauxie, are difficult to explain from the molecular point of view. 

Fig. 9a shows that in E. coli HfrSOOO, glucose is consumed before galactose. 
However , the mutant s train MM6, which contains a lesion in the PTS en- 



zyme I (jTanaka, et a 1.1 . 119671 ). consumes galactose before glucose (Fig. 9b). 
Likewise, E. coli strain 159 consumes mannitol before sorbitol (Fig. 9c), but 
the corresponding mutant strain 157, which contains a lesion in the PTS en- 
zyme II™*\ consumes sorbitol before mannitol (Fig. 9d). These phenotypic 
changes fall within the scope of the minimal model. In both mutants, the 
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Figure 9. Reversal of the diauxie in mutants of E coli: Upper pane l: (a) In strain 
HfrSOOO, glucose is consumed before galactose I Joseph et alJ . 1l98lh . The optical 



density (OD) shows a pronounced diauxie lag, regardless of the presence of cAMP 
in the culture, (b) In the correspond ing PTS-deficient strain, galactose is consumed 



before glucose (|Asensio et all Il963l ). Note that the evolution of the OD during the 



first 8 h is the same during growth on galactose [OD (GAL)] and glucose + galactose 
[OD (GLU+GAL)]. Furthermore, there is no consumption of glucose during this 
period. Lower panel: (c) In strain 159, mannitol (MTL) is consumed before sorbitol 
(SBL). (d) In th e corresponding enzvme II™*'-deficient strain, sorbitol is consumed 
before mannitol ( Lengeler and Liiil . 1972h . 



transport enzyme for the preferred substrate is impaired, i.e., Vs,2 decreases, so 
that Ki remains unchanged, but K2 increases and a decreases. If the changes in 
K,2 and a are sufficiently large. Fig. 6b implies that the substrate consumption 
pattern will shift from preferential consumption of 5*2 to preferential consump- 
tion of 5*1. 



It should be emphasized that the "reversal of the diauxie" is a natural conse- 
quence of the minimal model. This is because each enzyme inhibits the other 
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enzyme due to dilution by growth, i.e., the inhibition is mutual or competitive. 
Consequently, suppressing the uptake (and hence, the growth) on one of the 
substrates automatically tilts the balance of power in favor of the other sub- 
strate. In contrast, the "reversal of the diauxie" is difficult to explain in terms 
of molecular mechanisms alone. This is because in all the molecular mecha- 
nisms, the inhibition is unilateral rather than mutual. In E. coli, for instance, 
there are numerous mechanisms that allow PTS sugars, such as glucose and 
mannitol, to inhibit the synthesis of the enzymes for non-PTS substrates. But 
there is no mechanism for non-PTS substrates to inhibit the synthesis of PTS 
enzymes. This difficulty did not escape the attention of Asensio et al, who ob- 
served the reversal of the glucose-galactose diauxie (Fig. 9, top panel). Faced 
with the "reversal of the diauxie," they were compelled to conclude that the 
"diauxie is, at least in part, due to competitive effects at the permease level." 



3.1.2 Dependence of substrate consumption pattern on relative growth rates 

In order to consider the relationship between the substrate consumption pat- 
tern and the ratio of the single-substrate maximum specific growth rates, 
define 



~,max 
_ ^9,2 



'^9,1 

where r^f"" denotes maximum specific growth rate during single-substrate 
growth on saturating concentrations of Si. Now, the model implies that 



so that 



a r-^ = r, ag[Ki,K2) = 



£2 


lEoi 


a 











I-Bio 


-Hi + 


62 




—K2 + 



Ki + 4 



It follows that 

(1) The surface of ag{Ki. K2) separates the parameter space into two distinct 
regions: Above the surface, p > 1, i.e., r"^^ > r™^"", and below the surface, 

p<l, i.e.,r-r<r-r- 

(2) The surface of ag{Ki,K2) lies between the surfaces of a*{Ki,K2) and 
0£*{k,i, K2), i.e., 

Q;*(«;i, K2) < ag{Ki, K2) < a*{Ki, K2) (24) 

for all K,2 > (see Appendix C). Thus, the graph of a^, denoted by the 
dashed brown line in Fig. 6, lies between the graphs of a^, (blue curve) 
and a* (red curve). 
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Given these results, we can recast the classification of the substrate consump- 
tion patterns in terms of p. To this end, define 



P*(/«2) = 7 T = 1^2- 



ag{Ki,K,2) 2 
a*(Ki,K2) 1 2 



Then, there is preferential consumption of 5*1 (resp., S2) precisely when p < 
P*{i^2) (resp., p > p*{ki)), and simultaneous consumption if and only if p*(/t2) < 
p < p*(ki). Thus, and p* define the limits of p at which there is simultane- 
ous consumption. It turns out that p*(k2) increases from to 1 as K2 goes from 
to cx), and p*{i^i) decreases from 00 to as ni goes from to cxd (Fig. 10). We 
are now ready to discuss the relationship between the substrate consumption 
patterns and the ratio of the single-substrate maximum specific growth rates. 

The Harder & Dijkhuizen correlation states that when growth is diauxic, the 
preferred substrate is the one that, by itself, supports a higher maximum 
specific growth rate (p. 7). The model predictions are consistent with this cor- 
relation. This is already evident from Fig. 6: a < i.e., p < 1 in the region, 
a < a^,, corresponding to preferential consumption of Si, and a > Ug, i.e., 
p > 1 in the region a > a* corresponding to preferential consumption of 82- 
The same property is also manifested in Fig. 10, e.g., in the region, p < p^{k2), 
corresponding to preferential consumption of 6*1, p < 1 because the graph of 
p*(/«2) is always below 1. The manifestation of the Harder-Dijkhuizen correla- 
tion in this minimal model suggests that is an intrinsic property of the induc- 
tion and dilution kinetics. It can be explained without invoking goal-oriented 
regulatory mechani sms, which form the basis of in odels based on optimality 
princ iples (fKompala et a1.[ll98(Tl:lMaha,deva,n et a,1.l.l2002;.Rama,krishna et a,1 



Il99fih . 

Current experimental evidence suggests that the existence of comparable single- 
substrate maximum specific growth rates is, perhaps, necessary but not suffi- 
cient for simultaneous consumption (p. 8). However, Fig. 10 shows that this 
condition (p ^ 1) is neither necessary nor sufficient for simultaneous con- 
sumption. It is not necessary because when ki,K2 <C 1, there is simultaneous 
consumption for almost all p. It is not sufficient for simultaneous consumption 
because when Ki,K2 3> 1, simultaneous consumption is virtually impossible 
— it cannot be obtained unless p lies in a vanishingly small neighborhood of 
1. These results can be understood in terms of the limiting cases discussed 
above. If ki, K2 are small, the enzymes are quasi-constitutive, and they resist 
extinction, regardless of the maximum specific growth rates. As ki and K2 in- 
crease, the enzymes become progressively more vulnerable to extinction, and 
in the limit of large ni, k,2, they cannot coexist. 
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Figure 10. Classification of substrate consumption patterns in terms of the ratio 
of the single-substrate maximum specific growth rates. The full and dashed lines 
show the graphs of /9*(At2) and respectively. 

We note finally that unlike the standard Lotka-Volterra model for competing 
species, there are no parameter values that yield bistable enzyme dynamics 
(compare Figs. 5 and 6). We show below that bistability becomes feasible 
when the induction kinetics are sigmoidal. 

3.2 Case 2 (ni = 2, n2 = 1) 

In this case, the scaled equations are 



The key results, which are shown in detail in Appendix D, are as follows 

(1) The trivial steady, Eqq, always exists, regardless of the parameter values. 
It is always unstable. 

(2) The semitrivial steady state, £"10, exists if and only if Ki < 1, in which 
case it is unique, and given by 



dei 

m 

de2 



£2 

a 

H2 + ^2 




(ei + ae2) ei 



(ei + 0:62) €2 




It is stable (as a node) if and only if ei|g^^ exceeds the ei-intercept of the 
nontrivial nullcline for €2, i.e.. 
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Figure 11. The bifurcation diagram for the case ni = 2, 712 = 1 at any fixed K2 > 0. 
In the phase portraits, the nullchnes for ei and €2 are represented by full and dashed 
lines, respectively; stable and unstable steady states are represented by full and 
open circles, respectively. The graphs of a*, a*, ki = 1, and Og are represented by 
blue, red, green, and dashed brown curves, respectively. In the hatched region, 5*2 
is the preferred substrate for all preculturing conditions, even though it supports a 
maximum specific growth rate lower than that on Si. 

(3) The semitrivial steady state, Eqi, always exists, and is given by 



ei = 0, 62 



(26) 



It is always stable (as a node). 
(4) Nontrivial steady states exist only if ki < 1. Under these conditions, there 
are at most two nontrivial steady states. There is a unique nontrivial 
steady state if and only if 

< a < a*(/ti, K2), 

and it is unstable whenever it exists. There are two nontrivial steady 
states if and only if 



< Ki < /3 



2 + Ki 



\2(l + K 



a*(Ki, K2) < a < a {ki, K2), 



where q;*(ki, K2) is the value of a at which the nontrivial nullclines for 
61 and 62 touch. One of these steady states is stable and the other is 
unstable. 

(5) The surface of a*(Ki,K2) lies below the surface of a*{K,i,K,2) for all < 
Ki < P and K2 > 0. 



The bifurcation diagram shown in Fig. 11 implies the following classification 
of the substrate utilization patterns. 
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Figure 12. The bifurcation diagram obtained when a is increased at ki = 0.5 and 
K2 = 1. The red and blue curves show the semitrivial steady states, E'lo, -Boi, 
respectively. The green and black curves show the two nontrivial states (^^ii). The 
curves are full (resp., dashed) if the corresponding steady state is stable (resp., full). 
The bifurcation points are represented by full circles. 

(1) If a < a*, E\Q and Eq\ are stable, i.e., there is preferential consumption 
of 5*1 or 5*2, depending on the initial conditions. 

(2) If < Ki < /3 and a* < a < a*, i?oi and E^x are stable, i.e., there is 
preferential consumption of 5*2 or simultaneous consumption of 5*1 and 
5*2, depending on the initial conditions. 

(3) If /3 < Ki < 1, a > a* or Ki > 1, there is preferential consumption of S'2, 
regardless of the initial conditions. 

The surfaces of a* and a* define the locus of transcritical and fold (saddle- 
node) bifurcations, respectively (Fig. 12). If a is increased at any fixed < 
Ki < /? and K2 > 0, the substrate consumption pattern changes at a = 
from bistable dynamics involving preferential consumption of S\ or S'2 to 
bistable dynamics involving preferential consumption of S'2 or simultaneous 
consumption. This transition occurs via a transcritical bifurcation. At a = a*, 
the substrate consumption pattern switches to preferential consumption of S'2 
via a fold bifurcation. 

Comparison of Fig. 11 with Fig. 6 shows that certain features are preserved. 
Specifically, preferential consumption of S\ is feasible only at low a, and simul- 
taneous consumption occurs only if a has intermediate values and ^1,^2 are 
not too large. However, a unique property emerges in Fig. 11, namely, bista- 
bility. This is due to the sigmoidal induction kinetics for which ensure that 
preferential consumption of S'2 is feasible at all parameter values. 

It is also worth examining the relationship between the classification predicted 
by the model and the empirical classification based on the single-substrate 



26 



maximum specific growth rates. In this case 



Gl. I ^ 1^ ^ Gig 



2 



Now, > a* because — /t2 + \J K2 + ^ < 2/ K2 (see Appendix C). Furthermore, 
Q!g is zero at ki = 1. Thus, the graph of ag lies above the graph of a* (dashed 
brown line in Fig. 11). This implies that a substrate can be consumed prefer- 
entially even if it supports a lower maximum specific growth rate. Indeed, if 
the parameters lie in the region, a < ag, then 5*2 supports a lower maximum 
specific growth rate than Si, and yet, cells precultured on S2 consume this 
substrate preferentially. If the parameter values lie in the hatched region of 
Fig. 11, 5*2 is the preferred substrate, regardless of the manner in which the 
cells are precultured. 



3.2.1 Evidence of bistable substrate consumption patterns 

The bistable dynamics predicted by Fig. 11 have been observed in experiments. 

The bistable dynamics in the region, a^ < a < a* , correspond to preferential 
consumption of S2 if the preculture is grown on 5*2, and simultaneous con- 
sumption if the preculture is grown on Si. Two examples of this substrate 
consumption pattern were described in the Introduction, namelv. growth of 
P. aeruginosa on glucose plus citrate ( Hamilton and Dawesl . 195ji) , , I96Q. IQGlT) 



and gr owth of E. coli K12 on a mixture of glucose and pyruvate (lNara.ng et al 



1997bf ). Fig. 13 shows another example of this substrate consumption patte 



rn. 



When Streptococcus mutans GS5 is grown on a mixture of glucose and lactose, 
glucose-precultured cells consume glucose before lactose (Fig. 13a), whereas 
lactose-precultured cells consume both glucose and lactose (Fig. 13b). 

The bistable dynamics in the region, a < a^, correspond to preferential con- 
sumption of 5*1 if the preculture is grown on Si, and preferential consump- 
tion of S2 if the preculture is grown on 5*2. Furthermore, the maximum spe- 
cific growth on S2 is lower than that on Si. There is evidence suggesting 
the existence of this substrate consumption pattern. Tsuchiya and coworkers 
studied th e growth of Salmonella ttiyhimur i um o n a mixture of glucose and 
melibiose ( Kuroda et al. . 19921 : Okada et al. . 1981 ). They found that the wild- 



type strain LT2 consumed glucose before melibiose. However, the PTS enzyme 
I mutant, SB1476, yielded the bistable substrate consumption pattern corre- 
sponding to the region, a < a^,. Cells precultured on glucose consumed glucose 
preferentially (Fig. 13c), and cells precultured on melibiose consumed meli- 
biose preferentially (Fig. 13d). Moreover, the maximum specific growth rate 
on glucose (0.24 h~^) is significantly lower than that on melibiose (0.41 h~^). 
It should be noted that these experiments were done in the presence of 5 mM 
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Figure 13. Bistability in mixed-substrate growth: Upper panel: Growth 
of Stre ptococcus mutans GS5 on a mixture of lactose (LAC) and glucose 
(GLU) (|Liberma,n and Bleiwei^TlflsJ ). (a) Glucose is consumed preferentially if the 
cells are precultured on glucose, (b) Glucose and lactose are consumed simultane- 
ously if the cells are precultured on lactose. Lower panel: Growth of Salmonella 
typhimu rium SB1476 on a m ixture on a mixture of glucose (GLU) and melibiose 
(MEL) llKuroda et al .l. llflfl2h . (c) Glucose-precultured ce Is consume glucose before 
melibiose. (d) Melibiose-precultured cells consume melibiose before glucose. In (b) 
and (d), the concentration of glucose increases at i 8 and t 6 h, respectively. 
It is believed that is due to expulsion of the glucose produced from intracellular 
hydrolysis of lactose and melibiose, respectively. 



cAMP in the culture. However, at least in the case of glucose-precultured cells, 
the same phenotype was observed even in the absence of cAMP. 
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4 Discussion 



We have shown that a minimal model accounting for only enzyme induction 
and dilution captures and explains all the substrate consumption patterns 
observed in the experimental literature. In what follows, we discuss the ro- 
bustness of the model, and its implications for the problem of size regulation 
in development. 



4.1 Robustness of the model 



Given the simplicity of the model, it is necessary ask whether the properties 
of the model will be preserved if additional metabolic details and regulatory 
mechanisms are incorporated in the model. Now, the defining property of the 
minimal model is that the enzymes follow competitive dynamics. We show be- 
low that this property is not a consequence of the particular kinetics assumed 
in the model. It is the outcome of two very general characteristics possessed 
by most systems of mixed-substrate growth. 

To see this, it is u seful to consider the gen eralized Lotka-Volterra model for 
competing species (iHirsch and Sma1el . ll974 Chap. 12). This model postulates 



that the competitive interactions between two species are captured by the 
relations 

^ = m„N,), /i(0,iV2) = /2(iVi,0) = 0and |^,|^<0. 
dt 0N2 oNi 

In other words, the essence of competitive interactions can be distilled into 
two properties: 

(a) The growth of a species is impossible in the absence of that species 
{dNi/dt = whenever Ni = 0). 

(b) Each species inhibits the growth of the other species {dfi/dN2, df2/dNi) < 
0. 

These properties, by themselves, imply the existence of all the dynamics asso- 
ciated with competitive interactions, namely, the absence of limit cycles, and 
the existence of extinction and coexistence steady states. 

Now, properties (a) and (b) will be manifested in most systems of mixed- 
substrate growth. Indeed, the evolution of the enzymes during the first expo- 
nential growth phase can be described by the relations 

dci 

= gi{ei, 62) = re,i(ei, 62) - rg{ei, 62)6^. 
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(a) (b) 

Figure 14. Disappe arance of the diaux ie when the enzymes are not necessary for 
their own synthesis (|lna,da et a,l. . IflPfil ). Synthesis of /3-galactosidase persists in the 



presence of glucose if (a) lad is mutated so that lac transcription becomes consti- 
tutive or (b) IPTG is present in the medium. In the first case, lac transcription is 
constitutive, i.e., it persists even in the absence of the inducer. In the second case, 
lac transcription is no longer dependent on the existence of the permease. 

If we assume that 

(1) each enzyme is necessary for its own synthesis, i.e., r^^i = whenever 

Ci = 0, 

(2) each enzyme has either no effect or inhibits the synthesis of the other 
enzyme, i.e., dre,i/de2, dre,2/dei < 0, 

(3) the specific growth rate is an increasing function of Ci and 62, i.e., dvg/dei, 
drg/de2 > 0, 

then the enzymes satisfy both the hypotheses of the generalized model for 
competing species: (a) There is no enzyme synthesis in the absence of the 
enzyme {dci/dt = whenever Cj = 0), and (b) each enzyme inhibits the 
synthesis of the other enzyme {dgi/de2, dg2/dei < 0). Consequently, they will 
display extinction and coexistence dynamics. 

It remains to consider the generality of assumptions 1-3. 

Assumption 1 will be satisfied whenever the substrates are transported by 
unique inducible enzymes. In these cases, the enzymes are required for the 
existence of the inducer (cj = ^ = 0), and the inducers are necessary for 
the synthesis of the enzymes (xj = ^ re,t = 0); hence, Cj = ^ re,j = 0. 
One can imagine two cases in which assumption 1 is violated. First, if an 
enzyme is constitutive, it is synthesized even in the absence of the inducer 
{xi = =^ Te^i = 0). Second, in the presence of a gratuitous inducer, such as 
IPTG, which can enter the cell even in the absence of lactose permease, the 
enzyme is not required for the existence of the inducer (ci = ^ Xi = 0). In 
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both cases, the "extinction" steady state ceases to exist, and the substrates will 
be consumed simultaneously. This is consistent with experiments (Fig. 14). 



Assumption 2 will be satisfied provided the enzymes do not activate each other. 
But all the known regulatory mechanisms invariably entail direct or indirect 
inhibition of one of the enzymes by the other enzyme. This includes inducer 
exclusion (dephosphorylated enzyme IF^'^ inhibits lac permease), and cAMP 
activation (dephosphorylation of IP''^ causes a reduction of cAMP levels, which 
in turn inhibits lac transcription). 

Assumption 3 will be satisfied if the yield of biomass on a substrate dur- 
ing single-substrate growth does not change markedly during mixed-substrate 
growth. In the model, the yields were assumed to be constant. This is cer- 
tainly true for conservative substrate s since Yj = 1. It i s also observed to hol d 
in many mixtures of carbon sources ( Egli et al. . Il982l : Narang et al. . 1997bl ). 



However, it is conceivable that there are systems in which the yields vary with 
the enzyme levels. In such cases, the specific growth rate will have the form, 
rg{ei,e2) = yi(ei, 62)^,461 + ¥2(61, 62) Vs ^2^2- At present, the data is not suffi- 
cient for determining the extent to which the yields with the enzyme levels. 

It is therefore clear that inclusion of various regulatory mechanisms will en- 
hance the mutual inhibition due from dilution. However, the qualitative be- 
havior will be preserved, since the enzymes will still follow Lotka-Volterra 
dynamics. Thus, the key property of the model, namely, competitive dynam- 
ics of the enzymes, is quite robust insofar as the perturbations with respect 
to regulatory mechanisms are concerned. 

The notion that diauxic growth is the outcome of competitive interactions 
between the enzymes is not new. It can be found in the earliest papers on 
diauxic growth. In 1947, Monod noted that (|Monodl . Il947l . p. 254) 



"it appears that the mechanisms involved in diauxic inhibition have the char- 
acter of competitive interactions between different specific enzyme-forming 
systems." 



He observed, furthermore, that the (|Monodl . ll947l p. 259) 



"existence of competitive interactions in the synthesis of different specific 
enzymes appears to be a fact of fundamental significance in enzymatic adap- 
tation, and one for which any conception of the phenomenon should be able 
to account." 

However, these conclusions were based on the kinetics of the enzyme levels 
during diauxic growth, and had no mechanistic basis. 



The above argument, first made in iNarang l (|l998ar l. shows the mechanistic 
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basis of the competitive interactions in a mathematically precise fashion. 



4.2 Implication of the model for development 



Diauxic growth has played a critical role in shaping models of patterning in 
development. The first link between genetics and development was established 
in the late 40's by appealing to the following argument dGilbertl [20021. During 
diauxic growth, cells possessing identical genes synthesize different proteins 
at distinct times (namely, the first and second exponential growth phases). 
By analogy, patterning in differentiation could be viewed as the svnthesis of 
differ ent proteins at distinct times and locations ( Monodl . 1947 : Spiegelmanl . 



Il948i l. From this standpoint, diauxic growth and developmental patterning can 
be viewed as "temporal" and "spatiotemporal" differentiation, respectively. 

The subsequent discovery of the molecular mechanisms involved in develop- 
mental patterning have confirmed the above hypothesis. It has been found that 
developmental patterns are generated by mechanisms simil ar in principle, but 



more complex in detail, than those involved diauxic growth (jPtashne and Gann . 
20021 . Chap. 3). 



Despite remarkable successes in developmental patterning, there are outstand- 
ing questions about size regulation, i.e., the mechanisms bv which pattern- 
ing is coupled to grow th (|Dav and T^awrencl I2OOOI: iHafen and Stocked ^ ' 



Serrano and OTarrell 119971 ). Examples of such questions include: What de- 



termines the size of organs and organisms, i.e., why does their growth cease at 
a certain time? And why is development scale-invariant,, i.e., why is the size 
of the organs is proportional to the size of the organism? 

The model presented here may be relevant to the problem of size regula- 
tion. It shows that the "temporal" differentiation in the diauxie is coupled 
to growth, and this coupling is mediated by the process of enzyme dilution. 
Inasmuch as the diauxie is a paradigm of the mechanisms controlling cellu- 
lar differentiation, a similar mechanism may lie at the heart of the coupling 
between developmental patterning and growth. Based on the minimal model, 
one can speculate, for instance, that organ growth ceases at a certain time 
because growth-promoting enzymes are driven to "extinction" at sufficiently 
high growth rates. 

The model also has implications for the problem of scale invariance. In many 
mathematical models of development, pattern formation occurs when a ho- 
mogeneous steady state of a reaction-diffusion system 

dc 

— =DV^c-r{c,p) (27) 
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becomes unstable due to the onset of a Turing instability ()Murravl . {H 
Here, c{x, t) denotes the vector of morphogen concentrations, D is the matrix 
of diffusivities, and r{c,p) is the reaction rate vector expressed as a func- 
tion of c and a vector of parameters, p. In general, the patterns predicted by 
these models are not scale-invariant. However, this problem can be resolved 
if the system is fed more information about its size (say, L). For instance, 
perfect scale invariance is obtained if the diffusivities or rate constants are 
proportional to L^, and plausible mecha nisms for such a d ependence have 
been proposed dOthmer and Patel. i1 980: .Tshihara and Kaneka. 2W^) ■ 



In growing systems, however, information regarding the growth rate is con- 
stantly fed to the mechanism driving pattern formation. Indeed, in the pres- 
ence of growth, Eq. (27) becomes 



dc „ 
at 



DV'^c-r{c,p) -cV 



(28) 



where v{x,t) is the velocity vector field, v ■ Vc is the accumulation of the 
morphogens due to convection, V ■ f is the specific growth rate, and cV ■ v is 
the dilution of the morphogens due to growth. Crampin et al have shown that 
these equations exhibit a certain degree of scale invariance — as the system 
grows, the number of pattern eleinents remain s the same despite a doubling 
of the system size ( Crampin et al. . 1999L 20021 ). Further analysis of this class 
of equations offers the promise of deeper insights into the coupling between 
patterning and growth. 



5 Conclusions 

(1) We showed that a minimal model accounting for enzyme induction and 
dilution, but not cAMP activation and inducer exclusion, captures and 
explains all the observed substrate consumption patterns, including di- 
auxic growth, simultaneous consumption, and bistable growth. This sug- 
gests that the dynamics characteristic of mixed-substrate growth are al- 
ready inherent in the minimal structure associated with induction and 
dilution. We find that many of the molecular mechanisms, such as in- 
ducer exclusion, serve to amplify these inherent dynamics. 

(2) We constructed bifurcation diagrams showing the parameter values at 
which the various substrate consumption patterns will be observed. The 
bifurcation diagrams explain the phenotypic responses to various genetic 
perturbations, including lesions in the genes for the repressor, operator, 
and the transport enzymes. Importantly, they provide a simple explana- 
tion for the "reversal of the diauxie," a phenomenon which is quite difficult 
to explain in terms of molecular mechanisms. The bifurcation diagrams 
also provide deep insights into the mechanisms underlying the empiri- 
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cally observed correlations between the substrate consumption patterns 
and the single-substrate growth rates. We found that 

(a) When the induction kinetics are hyperbolic, the preferred substrate 
is always the one that that supports a higher growth rate. This cor- 
relation is, therefore, unlikely to be the outcome of optimal design. 
It is a natural consequence of the fact that the enzymatic dynamics 
are governed by the rates of induction and dilution. 

If induction is sigmoidal, it is possible for the preferred substrate to 
support a lower growth rate than the less preferred substrate. We 
presented experimental data illustrating this case. 

(b) The existence of comparable growth rates is neither necessary nor 
sufficient for simultaneous consumption. When the saturation con- 
stants are small, simultaneous consumption occurs regardless of the 
maximum specific growth rates, since induction is quasi-constitutive. 
If the saturation constants are large, simultaneous consumption is 
impossible even if the growth rates are comparable. 

(3) The key property of the model, namely, competitive dynamics of the 
enzymes, is quite robust with respect to structural perturbations. 

(4) The model may have implications for the problem of size regulation in 
development, since it provides a mechanism for coupling differentiation 
and growth, namely, protein dilution. 



Acknowledgements 



This research was supported in part with funds from the National Science 
Foundation under contract NSF DMS-0517954. 



A Interpretation of rii, j, and Ke^i in terms of molecular interac- 
tions 



To express n^, K j, and K ^ j in terms of rnolecul ar parameters, we appeal to 
the Yagil & Yagil model ( Yagil and Yagil 1971 ). For notational clarity, we 
shall ignore the subscript i for the substrate; thus, the operator, inducer, and 
repressor will be denoted by O, X, and R, respectively. Furthermore, their 
concentrations will be denoted [O], [X], and [R], respectively. 



The Yagil & Yagil model views induction as the outcome of a competition for 
the repressor between the operator and the inducer. Induction occurs when 
the inducer molecules sequester the repressors away from the operator. The 
competitive interactions are represented by two binding equilibria 
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R + O^R-0, Ko 
R + nX ^R- Xn 

where n denotes the number of inducer molecules that bind to 1 molecule of 
repressor; [R ■ O], [R ■ X„] denote the concentrations of the complexes, R ■ O, 
R-Xn, respectively; and Kg, denote the dissociation constants for the two 
equilibria. 



[R-oy 

\R-XS 



(A.l) 
(A.2) 



It is assumed that 



(1) Enzyme synthesis is limited by the transcription rate, i.e., translation 
is not limiting. Thus, the specific enzyme synthesis rate is proportional 
to the specific transcription rate. Furthermore, the specific transcription 
rate is proportional to the concentration of the free operator, i.e., 

r.^vpi (A.3) 

where v denotes the enzyme synthesis rate per unit mass of operator. 

(2) The total concentrations of O and R, denoted [0]t and [R\t, respectively, 
are conserved, i.e., 

[0]t=[0] + [R-Ol (A.4) 
[R]t^[R] + [R-0] + [R-Xn]. (A.5) 

These two relations, together with Eqs. (A.1)-(A.2), constitute 4 equa- 
tions in 4 unknowns, namely, [O], [i?], [i?-0], and [R • In principle, 
these equations can be solved for [O], and substituted in (A.3) to obtain 
Tg. However, since the solution is cumbersome, it is convenient to make 
the following additional assumption. 

(3) The repressor is bound primarily to the inducer (rather than the opera- 
tor), i.e., 

< [R-Xn]. 

This assumption is valid under most conditions because the operator 
concentration (~2 per cell) is significantly smaller than the inducer con- 
centration. 



These assumptions yield 

""^ " ""^^ ~ ^'K^{l + [R]t/Kg) + [Xr 

where = i^[Ot\. 

In the case of constitutive enzymes, the repressor has a weak affinity for the 
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operator, i.e., Ko ^ [-R]t, so that 



Te ^ Ve 

regardless of the inducer concentration. 

In the case of inducible enzymes, the repressor has a high affinity for the 
operator, i.e., Ko -C [R]t, so that 

+ [xr 

Te ~ 



■K4R]t/Ko+[X]^' 

This is a monotonically increasing function of [X] with a small nonzero inter- 
cept. Neglecting this small "basal" enzyme synthesis rate yields^ 



B Derivation of the equations 



Equations (3-5) implicitly define the specific growth rate and the evolution 
of the cell density. To see this, observe that since all the intracellular con- 
centrations are expressed as mass fractions (g/gdw), their sum equals 1, i.e., 
xi + X2 + ei + 62 + c~ = 1. Hence, addition of equations (3-5) yields 

^ , Idc 

i=i i=i ^ 

which can be rewritten in the more familiar form 
dc ^ ^ 

i=l i=l 

where Vg denotes the specific growth rate. 



We c an simplify the model by observing that Xi ~ 10 ^ g/gdw ( Chung and Stephanopoulo^ . 



Il99fih and rs,i,rx,i ~ 1 g gdw~^ h~^. Thus, Xi attains quasisteady state on a 
time scale of 10~^ h. Moreover, the dilution term VgXi ~ 10~^ g gdw~^ h^^ 
is negligibly small compared to Tg^i^rx^i. Hence, within a few seconds, (3) be- 
comes, ~ Ts^i — r^^i, which implies that r„ ^ J2i'Yi^s,i, i.e., Yi is essentially 



^ Recent evidence suggests that in the case of the lac operon, the cooperativity does 
not arise from the binding of two inducer molecules to a single repressor molecule. 
Instead, it mig ht be due to the coo perative binding of a single repressor molecule to 
two operators (jOehler et al. 
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the yield of biomass on Si. Thus, we arrive at the equations 



dc 

dsj 

dt 

dcj 

dt 

Xi 



re,i - {yirs,i + y2rs,2)ei, 

Vs,ieiSi/{Ks,i + Sj) 
k ■ 



C~ ^ 1 — Xi — X2 — ei — 62 



where Xi is obtained by solving the quasisteady state relation, r^i^ Vg 



C Stability analysis of case 1 (ni = ^2 = 1) 



In this case, the steady states satisfy the equations 







1 



— ei — ae2 ei, 



= ( a ei — ae2 ] £2. 

K2 + £2 



and the Jacobian at any (ei, €2) is 



-^(ei,e2) 



(Ki+ei)2 



— 2ei — ae2 



—aei 



-62 



(k2 



^-e,-2ae2 



C.l Trivial steady state 



It is evident that Eqq exists, regardless of the parameter values. It is always 
an unstable node since 

' ^ 0" 
^ 

which implies that both eigenvalues, l/ni and 1/^2, are positive. 
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C.2 Semitrivial steady states 



The semitrivial steady state, £^io, always exists. It is unique and given by 



-, €2 = 0. 



Since 1/(aci + ei) = ei at £'10, the Jacobian at this steady state is 



Kiei 



2t\ — aei 



--ex 



and the eigenvalues are 

Ai = zei = — ei < U, X^ — e\. 

Ki + ei Ki + ei 

Hence, £"10 is stable (as a node) if and only if 



\i = ei<0 <S^Q;<ft;2eiF ^ — • 

K2 a K2 



(C.l) 



A similar analysis of the semitrivial steady state, £01, shows that it always 
exists. It is unique and given by 



ei = 0, £2 

It is stable (as a node) if and only if 

1 



-K2 + v^^iTI 



OL > 



(C.2) 



C.3 Nontrivial steady state (s) 



The nontrivial steady state(s), En, satisfy the equations 



= 



1 



— ei — a€2. 



1 

= a ei — ae2- 

K2 + €2 

Eliminating £2 from these equations yields the equation 

1 



Ki + ei 



[1 + a )ei + a(aK,i — K2), 



(C.3) 
(C.4) 
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which has at most 1 positive root, and it exists if and only if 
a{aKi — K2) < — Eqi is unstable. 

On the other hand, eliminating ei from (C.3-C.4) yields the equation 

a ,1 . K2 

= ( — h 0()e2 H Ki 



K2 + €2 a a 
which has at most 1 positive root, and it exists if and only if 

Ki < — -v^ hio is unstable. 

a K2 

Thus, Ell exists if and only if both £^10 and Eqi are unstable. Furthermore, it 
is unique whenever it exists, and is given by 



{aK2 — Ki) — 2a'^Ki + {aK2 — + 4 [(1 + a^) + aK2 
" 2(1 + q;2) 

-a{aK2 - Ki) - 2k2 + {aK2 - kiY + 4 [(1 + a^) + aK2 



62 = 



2(1 + q;2) 

It turns out that En is stable whenever it exists since (C.3-C.4) imply that 

-^11 = 7 — xn - 2ei - ae2 = — (^i + 2ei + aea) < 

[Ki + eiY Ki + ei 

and 

J22 = 7 ^ - ei - 2ae2 = {aK2 + ei + 2ae2) < 

(^2 + 62) «2 + ^2 

SO that tr J {En) < and 

det J {Ell) — + OiK'2) €1 + a {ki + aK2) €2 + ^aeie2 + 2Q;^e2 > 0. 
Hence, the eigenvalues of J{Eii) have negative real parts. 
We conclude that En exists if and only if 

Both £'10 and £'01 are unstable <S=> j — < a < K2 ei\^ . (C5) 

Furthermore, it is unique and stable whenever it exists. 
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C.4 Disposition of the surfaces of oc^, a*, and ag 

The surface of ag lies between the surfaces of q;*(«;i, K2) and a*{K\, K2), i.e., 

Oi*{K,i, K2) < ag{Ki, K2) < a*{Ki, K2) (C.6) 
for all Ki,K2 > 0. To see this, observe that 

/I x\'^ x'^ -X + \/x'^ + 4 1 

U + ^' + T^ 2 <J 

for all a; > 0. Hence 

-Ki + y'ftf + 4 -Ki + ^k\ + 4: ^ 2/ki 

for all Ki,K2 > 0, and (C.6) follows immediately from the definitions of a*, 
a*, and ag. 



D Stability analysis of case 2 (ni = 2, n2 — 1) 



In this case, the steady states satisfy the equations 



= f 2^] 2 - ^1 - '^^2 ) ei, 



= a- 



K2 + ^2 

and the Jacobian at any (ei, 62) is 



- ei - ae2 62- 



J{ei,e2) 



2Kfei r, 

R+i^ - 2ei - aes 



-62 



aK2 



(k2+€2)2 



-aei 



ei — 2ae2 



In what follows, we study the conditions on the parameter values for the 
existence and stability of all four types of steady states. 



D. 1 Trivial steady state 



The trivial steady, £'005 always exists, regardless of the parameter values. The 
Jacobian is singular at this steady state, but we can infer its stability from the 
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dynamics on the invariant lines, ei = and 62 = 0. Indeed, in the neighborhood 
of -Boo, 



dei 
Hi 



..0 'VI J' dt 



a 



ei=0 



— €2 > 0. 

/«2 



Hence, -Eqo is a nonhyperbolic saddle if < 1 and a nonhyperbolic unstable 
node if Ki > 1. 



D.2 Semitrivial steady states 



The semitrivial steady state, £'10, exists provided < 1, in which case it is 
unique, and given by 



ei = ^1 - Kf, 62 = 0. 
Since el + k^i = 1 at this steady state, the Jacobian is 



2 {kI — 1) ei —aei 







K2 



- ei 



and the eigenvalues are 



Ai = 2 (kI - 1) ei < 0, A2 = ei. 

Hence, £^10 is stable (as a node) if and only if 



cx I 

A2 = ciIb <Q < K,2yl - K?. 

K2 



(D.l) 



Analysis of the steady state, £01, shows that this steady state always exists, 
and is given by 



-K2 + V«2 + 4 

ei = 0, 62 = ^ 

The Jacobian at this steady state is 

—ae2 

-62 ^^-2a62 

and the eigenvalues are 

\ / n \ '^2 + 262 

Ai = — q;62 < U, A2 = — Q;e2 < U. 

K2 + 62 

We conclude that £01 always exists and is stable (as a node). 



(D.2) 
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D.3 Nontrivial steady state (s) 



D. 3. 1 Existence 

The nontrivial steady states satisfy the equations 

= o ^ o — ^i — cte2, (D-3) 
= a ei — ae2. (D.4) 

If Ki > 0, there are no nontrivial steady states, since (D.3) cannot be satisfied 
for any ei, £2 > 0. Indeed, 

ei / 1 \ 

Ki> 1 ^ —2 2 ~ "^^2 = ^1 —2 2 ~ ■'- ~ "^^2 < 

for all ei, 62 > 0. Henceforth, we shall assume that < < 1 and K.2 > 0. 

We begin by introducing a change of coordinates that reduces the problem 
to the existence of roots on a finite interval. Letting 62 = I'ei, we rewrite the 
above system as 

1 

= ei(l + av). 



K2 + ve-i 

Solving the first equation for ei, and rewriting the second equation, we obtain 



ei = \l T-r '^i' 

a 2 



1 + va 

Hence, we obtain the following equation for u: 



a f 1 o\ /I 9 
- + i^2\ — «1- 



l-\-va av J \ 1 -\- av 

Multiplying through by (1 + av), we obtain 



a — v 



1 - + av) \ + «;2\/(l + ai/) [1 - + av)\. 



Finally, we let 

z — K\il -\- av) , V — — ( — 1 

a \ Ki 



42 



and rewrite the above equation as 

1^-^{z-kI){1-z) + ^^z{1-z) = f{z, a). (D.5) 

Observe that z uniquely determines u, and hence both ei and €2. To ensure 
that both Si > 0, z must belong to the interval (k^, 1). Thus, the problem is 
reduced to the existence of roots of (D.5) on the finite interval {nf, 1). 

Several properties of / are immediate: 

/(l,a) = 0<l, f{Kla)^^ 

a 

and 

, , , 2 [4z(l - + (1 - 2^)2] 

so that f{z, a) is strictly concave down in z for each a > and for all < 
z < 1. In particular, this implies that (D.5) has a unique root z e [k^, 1] if 
f{K,l,a) > 1, and at most two roots ii f{Kl,a) < 1. Furthermore, we observe 
that fz{z, a) — > —00 as z — > 1, 

/^K, a) = , , 

2{1 - kI) + 4(1 - 2kI) 



and 



Hence, /^(ki, a*) > if and only if 2(1 — nf) + kI{1 — 2Kf) > 0, i.e.. 



Ki < P = 



Kin 



\2{l 



4) 



Finally, we observe that f{z, a) decreases in a for each fixed z, and /(z, a) — > 
as a — > +00. 

We conclude that there are two scenarios 

(1) If Ki < (3, then there exists a* > a* such that (D.5) has only one root in 
(k^, 1) for all a < a*, two distinct roots in (k^, 1) for all a* < a < a*, 
and no roots for a > a*. 

When a — a^:, (D.o) has two roots: z — and another root in {nf, 1). 
When a = a*, (D.5) admits a double root in {nf, 1). 

(2) If Ki > P, then (D.5) admits only one root in {kI, 1) for all a < a*, and 
no roots for a > a*. 

A single root z — k\ occurs if and only if a = a*. 
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Thus, we obtain the bifurcation diagram shown in Figure 11. 



D.3.2 Computation of a* {hii, K2) 



We have shown above that a* exists for all ki < /3. Here, we present an 
algorithm for computing a*, which rests upon the fact that the nontrivial 
nuUclines for ei and 62 touch at a = a*. 

Eliminating €2 from (D.3-D.4) yields the quartic polynomial 



■ (1 + a^) - 4 ■ otH2 + ■ 2a^K\ - (1 



0. (D.6) 



This equation has two equal real roots if and only if the discriminant is 
zero (IDicksonl . Il914i p. 41), i.e.. 



A 



a 



Co (a^) + Cl {a^^ + C2 (a^) + C3 







(D.7) 



where 

Co 

Cl 
C2 

C3 



KlU + Ki) >0, 



16 



(k^ - 4) - {kI + 2) (t^l + 4) - 4/t? (k^ + 4) (k^ + 8 
16/t^ + (k^ - e) + (k^ + 4)^ - 4k2 + 4) (bKl + 12 



iik^ -44 



f (i-K;n4(i-.;) + 



> 0. 



For every Ki, K2 > 0, equation (D.7) has three nonzero roots. One of these 
roots is negative since cq and C3 are positive. Computations show that the the 
remaining two roots are also positive. However, the nullclines touch in the first 
quadrant, €1,62 > 0, only if ki < /?, and a is the largest positive root. Thus, 
a*{K,i, H2) is the largest of the three roots of (D.7). 



D.3.3 Stability 

The stability of the steady states follows from the geometry of the nontrivial 
nullclines for ei and 62- Indeed, it is known from the theory of the generalized 
Lotka-Volterra model for competing species that a nontrivial steady state is 
stable if and only if in the neighborhood of the nontrivial steady state, the 
nontrivial nullclines for both ei and €2 have negative slopes, but the slope of the 
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nontrivial nullcline for €2 is more negative (i.e., higher in absolute value) than 
the slope of the nontrivial nullcline for ei (jHirsch and STna1el . ll974 Chap. 12). 
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